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Improved  Imaging 
of  Space  Objects 


Pristine 


Good  Seeing  Marginal  Seeing 

Terminator  Terminator 


Daytime 


Simulations  of  the  Hubble  Space  Telescope  as  it  would  appear 
from  the  3.6  m  AEOS  telescope  at  a  range  of  700  km  in  1  ms 
exposures  at  0.9  pm  wavelength  under  a  range  of  seeing 
conditions. 
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V  4  Improved  Extraction  of  Information 

from  ground-based  SSA  imagery 


Current  capability  New  capability 


Current  capability  using 
data  from  3.6m  alone 
(D/r0=27) 


Proposed  capability 
using  data  from  both 
the  3.6m  and  1.6m 


Leveraging  AF  telescopes  in  proximity 

•  Add  Resolution  diversity  (multiple  telescopes) 

•  Enforce  intra-  and  inter-channel  consistency 

•  Maintain  high-resolution  under  strong 
turbulence 


Potential  to  push  SSA  capabilities  to 
severe  seeing  conditions  (D/r0>  30) 


Stuart  Jefferies,  U  Hawaii 
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W  Augmen,ed  Hybrid  Input-Output 

Phase  Retrieval 


Truth  Image,  g 


Simulated  LR 


True  Image 


f 

g^G 


Interferometer  Array 


Imaging  Telescope 


ll/IT  Reconstruction 


Phase  Retrieval 
Algorithm 


Object  Estimate 


g 


Peter  Crabtree,  AFRL/RVB 
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Local  Maxima  Structure  of 
Wavefront  Estimation 


Each  speckle  associated  with 
subset  of  pupil. 

Each  maximum  corresponds  to  one 
mapping  of  subsets  to  speckles. 


Provides  physical  description  of  the  local  maxima  associated  with  speckle 
imaging  estimation  problems 


Allows  properties  of  local  maxima  to  be  derived  from  the  Kolmogorov 
model  of  atmospheric  turbulence: 


Expected  distance  between 
local  maxima 

0.2 
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0.1 

0.05 

5  10  15  20  25  30 
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Chuck  Matson,  AFRL/RDS 
Brandoch  Calef,  Boeing 


Optimal  number  of  modes 
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Wavefront  estimation  error 
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V  4  Algorithms  for  Multi-Frame  Blind 

Deconvolution 


Initial  PSF  Estimates  using  Wave  Front  Sensor 


■Properly  defined  constraints  to  define 
objective  functions 


Truth  New  FF  alg  Standard  alg 


Accelerating  MFBD  Convergence 


■Partial  second  derivative  preconditioning 

■  Implicit  Filtering  (C.T.  Kelley,  SIAM,  2011) 
to  avoid  local  minimum  trap 


o.i 


0.05  \ 

°  10  20  30  40  50  60  70  80  90  100 

_  iteration 

James  Nagy,  Emory  University 
Stuart  Jefferies,  University  of  Hawaii 


*  Faster  convergence  with  higher 
likelihood  of  achieving  the  global 
minimum. 

*  Improved  quality  of  image  restorations. 

*  Extended  limits  of  SSA  imaging 
capabilities. 
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Challenge: 

-Limited  information  on  each  object 

-Few  sensors  for  a  big  space 

-Cost  of  big  telescopes,  space-based  sensors 

-Lack  of  a  priori  knowledge 

-Potentially  massive  data  sets 


Questions: 

-How  to  extract  information  from  an  unresolved  object 
-How  to  search  large/sparse  data  sets 
-What  information  is  needed,  what  is  not  needed 
-Are  there  un-tapped  sources  of  information 


Path  Forward: 

-Use  of  smaller,  cheaper,  more  diverse  sensors 
-Sensor/Information  fusion 
-Information  theory 
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Satellite  Superquadrics 
Shape  Estimation 


ei,ea,e3  >  0 


* 

^  0  *4 

4  • 


Different  super-ellipsoids 
from  changing  sl5  s2,  s3 


•Estimate  rotational  / 
translational  motions 
from  data  via  Fourier 
Descriptors 

•Combine  with  body 
silhouettes,  estimate 
body’s  3D  convex  hull 

•Estimate  superquadric 
parameters. 


0* 


Superquadrics  from  computer  vision 
simulate  complicated  surfaces 


Sudhakar  Prasad,  UNM 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


9 


w 

<«» 


S/C  Identification  from  Surface 

Characterization 


Attitude  Model 


Wire-frame  Shape  Model 


Paul  Kervin,  AFRL/RDS 


Candidate  Material  BRDFs 


Whole-body,  Multi-band 
Observations 


\  <1 


11-band  OCS 
No  noise 


* 


Object  =  22176  Site  =  AEOS 


11-band  OCS 
10%  Gaussian  noise 


12  14  16  18  20  22  24  26  28 

Minutes  past  2009-01-02  15:00:00 


Best  Material  Identifications 
for  Satellite  Components 


1  2  3  4  5  6  7 

Material  Number 
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What  is  minimum  data  set  for 
satellite  identification? 


Attitude  Model 


Wire-frame  Shape  Model 


Whole-body,  Multi-band 
Observations 


No  material 
BRDF  library 
employed  in 
retrieval 


11-band  OCS 
No  noise 


* 


Object  =  22176  Site  =  AEOS 


11-band  OCS 
10%  Gaussian  noise 


12  14  16  18  20  22  24  26  28 

Minutes  past  2009-01-02  15:00:00 


Retrieved  BRDFs  for  all 
detected  satellite  components 


Log-Scale  BRDF  Surface  for  Component  1  (K  =  0.6  p.m) 
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Coherent  Sensing  in  the  Presence  of 
Diffusing  Backgrounds  or  Obscurants 


target 


P  -  uniformly 
polarized  but 
usually  speckled 
(spatially  partially 
coherent) 


Field  available  to 
sensor  is  the 
coherent 
superposition  of 
both  contributions 


Step  1 


Extract  the  state  of  polarization  by 
mixing  the  fluctuating  field  with  an 
uncorrelated  reference  having  the 
same  mean  optical  frequency. 


Step  2 


Determine  the  spatial 
distribution  of 
polarization  states  of 
the  measured  field 


Step  3 


Map  of  similarity  of  polarization  states 
wrt  a  reference  (Complex  Degree  of 
Mutual  Polarization) 


An  analogous  situation  exists  for  spatially  unresolved 
signals  that  vary  in  time.  A  similar  sensing  procedure  can 
be  developed  where  the  spatial  maps  are  replaced  by 
time  series  of  polarization  resolved  measurements 


*  Relative  strength  of  field  lP  can  be  found  from  the  global 
degree  of  polarization;  this  is  the  target  reflectivity 

*  Spatial  properties  of  lP  relates  to  target  morphology  - 
size,  shape,  surface 


Aristide  Dogariu,  UCF 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


12 


Accurately  Predicting  the 
Location  of  Space  Objects 


M&nihly  N  umber  of  Objects  in  Earth  Orbit  by  Object  Type 


Cornerstone  capability  to  all  SSA 


Challenges: 

-  Large  number  of  objects  (>  20,000  of  size 
greater  than  10cm) 

-  Limited  observations  from  multiple 
sensors 

-  Non-Gaussian  errors,  uncertainty 
analysis 

-  High  fidelity  models  representing  space 
weather 

-  Conjunction  analysis  (combinatorial 
problem) 

-  Closely  spaced  objects 

-  Data  association  from  diverse  sensors 

-  Optimal  sensing  for  maximum 
information  gathering  and  hazard 
assessment  in  a  timely  manner 

-  Efficient  computations 
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Accuracy  of  Orbit  Predictions 


Current  uncertainty  assumptions  are 
Gaussian  (e.g.  constrained  to  error 
ellipses) 

•  AFSPC  has  shown  that  these  uncertainties 
are  often  unrealistic 

•  Uncertainty  has  wrong  size  and  shape 

Non-Gaussian  uncertainties  exist 

•  Objects  that  are  tracked  infrequently  or 
detected  for  the  first  time 

SSA  requires  realistic,  quantifiable,  and 
usable  measures  of  orbit  uncertainty 
(ambiguity) 

•  Correct  knowledge  leads  to  meaningful  and 
appropriate  UCT  mitigation,  collision 
probability  computations,  change  detection, 
etc. 


2.797  2.798  2.799  2.8  2.801 

^-Position  [km]  xl° 
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Accuracy  of  Orbit  Predictions 


Adaptive  Entropy  Gaussian  Information 
Synthesis  (AEGIS) 

•  Gaussian  Sum  approximations 

•  Information-Theoretic  measures  of  ambiguity 

Dr.  MoribaJah  AFRL/RV 


-6981  -6980  -6979  -6978  -6977  -6976  -6975 
x(in  km) 


Generalized  Polynomial  Chaos  (gPC)  +  Adaptive 
Gaussian  Mixture  Model  (AGMM) 

•  Pose  the  optimal  information  collection  problem  as  a 
stochastic  control  problem 

•  Performance  and  robustness  metrics  are  derived  in  from 
information  theory. 

Dr.  Puneet  Singla,  U  Buffalo  (YIP) 

/vir- 


Multiple  Hypothesis  Tracking  (MHT) 

•Nonlinear  filter 

•Complexity  of  the  Unscented  Kalman  Filter  with  the  performance  of  a 
Gaussian  Sum  Filter 

Dr  Aubrey  Poore,  Numerica  Corp. 
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Rigorous  Characterization  of 
Thrusting  Spacecraft 


(XfPf) 


Optimal  Connecting  x*(t) 
with  maneuvers 


‘ _ Propagated  x(t) ,  P(t) 


Control  Distance: 


d 


=  /|u*(r)|  dr  ^ 


Orbit  Range  in  semi-major  axis,  eccentricity 
and  inclination  for  a  satellite  starting  in  GTO 


Objective 

•A  rigorous  and  hypothesis-free  approach  to  characterizing 
the  actions  of  a  thrusting  satellite 

•Independent  of  whether  it  follows  an  impulsive,  low-thrust, 
or  natural  force  augmentation  approach. 


Approach 

•Application  of  Optimal  Control  Theory 
•Determine  the  limiting  control  expenditures  to  yield  a 
change  in  state 

•Discriminate  between  natural  forces  and  artificial  controls 
Delimit  “orbit  range”  accessibility  based  on  fuel  costs. 

DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


Surface  encloses  all  possible  orbit  states 
reachable  with  a  fuel  budget  of  1  km/s. 


D.J.  Scheeres  (University  of  Colorado) 
K.T.  Alfriend  (T exas  A&M  University) 
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Low-luminosity  Objects 


-192  -128  -  64  0  64 


X  (pixel) 


Yukihito  Kitazawa,  IHI  Corp 
Toshiya  Hanada,  Kyushu  U 


•  Object  with  luminosity 
lower  than  the  limiting 
magnitude  of  a  single 
image  is  successfully 
detected  with  a  set  of  32 
successive  images 

•  Combination  of  the 
Stacking  Method  with 
Motion  Prediction  reduces 
computation  time  by 

a  factor  of  35 
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Propagation  through  Deep 
Optical  Turbulence 


DE  systems  can’t  function  across 
all  parts  of  the  battlespace. 


Current  efforts: 

MURI 
JTO  MRI 
STTR 
YIP 

Conventional  grants  (4) 

Lab  Tasks  (2) 

AFIT  (2) 

International  research  grant 
NATO  Study 
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Horizontal 

Characterization 


Propagation 
and  Compensatio 


Propagation  Range  Green  Beacon  on  AEOS  Laser  beam  from  Haleakala 


Experiment  Description 

•  3  laser  beams  at  0.532,  1.06  and  1.55  microns 

•  Mauna  Loa  on  Hawaii  (11,100  ft)  to  Haleakala  on  Maui  (10,500  ft) 

•  Path  range  is  149  km 


Observed  Deviations 

•Lack  of  stationarity  and  isotropy  in  the  data 

•  Random  motion  of  speckles  and  speckle  pattern 

-violates  Frozen  flow  hypothesis 

•  Presence  of  unusually  large  amplitude  spikes 

-(coherent  structures) 

•  3  to  5  times  larger  measured  intensity  variance 

•  Aperture  averaging  is  not  effective  for  large  apertures 

•  Failure  to  scale  properly  with  wavelength  changes 


•  Approach  1:  Wave  propagation  through  conventional  Anisotropic  non-Kolmogorov  turbulence 

•  Good  match  over  100  km  paths  at  moderate  turbulence  between  simulations  and  analysis 

•  Phase  screen  simulations  do  not  predict  the  COMBAT  data  behavior 

•  Anisotropy  of  the  atmosphere  is  of  some  influence  but  is  not  a  major  factor. 


•  Approach  2:  Use  the  Nosov  turbulence  spectrum  with  phase  screen  simulations 

•  Spectrum  includes  possible  incipient  turbulence  and  coherent  structure  formation  in  open  air 

•  Non-K  exponents  unrestricted,  inner  scales  >  1  cm  and  outer  scales  <  1  meter 

•  The  substantial  difference  in  the  correlation  functions  at  larger  spacing  not  expected 


Rao  Gudimetla,  AFRL/RDS 
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Branch  Points  indicate  the  presence  of  Orbital 
Angular  Momentum  (OAM) 


OAM  pairs  can  be  experimentally 
created  on-demand 


DM  commands 

(hidden  phase) 

• 

J 

Deep  Turbulence  (branch  points)  can  be  expunged 


disturbance 


(hidden  phase) 


atmospheric  branch  points 

atmospheric  branch  points 

final  result 


Headway  into  a  problem  previously  thought  to  be  intractable 


Darryl  Sanchez,  AFRL/RDS 
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Propagation  of  Polarization  Modulated 
Beams  through  a  Turbulent  Atmospher 


SIT 


SUSCEPTIBILITY  TENSOR 
*1 


Absorption 

G 


Stimulated  Emission 


P  =  6 


0 


EIT 


11 

A 12 

X,  13 

21 

*  22 

*23 

31 

*32 

*33 
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Probe  beam  frequency 

SIT  &  EIT  well  known,  but  pulse 
(transient)  effects 


Tensor  nature  of  susceptibility/refractive  index/absorption/refraction: 
Describing  anisotropic  media 


sigl  orthogonally  polarized  to  (sig2  +  sig3) 


Poincare  Sphere  Representation 


Polarization  modulation  (POLMOD)  for  CW 
beams:  Transient  in  polarization  compatibility, 
but ,  but  using  CW:  POLMOD  achieves  CW-SIT 
&  CW-EIT 


POLMOD 
achieved  by 
combining  two 
orthogonally 

polarized  beams  of 
time:  1000  samples  different 

=  1  psec. 

wavelength 


Terrance  Barrett,  BSEI 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


Transitions 


•STTR  -  Gravity  model  reformulation 

•Accurate  high-speed  gravity  model  for  AFSPC 
•STTR  -  Prediction  of  Satellite  Ballistic  Coefficients 
•Increased  accuracy  of  drag  prediction  for  AFSPC 
•STTR  -  Synthetic  scenery  generation 
•Includes  obscurants 

•STTR  -  Beam  control  for  optical  phased  arrays 

•Increased  beam  quality  and  higher  power,  extended 
beacons 

•JTO  MRI  -  Airborne  Aero-Optics  Laboratory 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  is  unlimited. 


22 


New  Data  from  the  AAOL 


As  of  July  2010 


All  from  Wind  Tunnel 


By  February  2011 
A  Sample  of  Flight  Data 


E 

o 


Q 

CL 

O 


a 


Morgan  et  al.,  2011 


CFD  confirmation  and  analysis 


Eric  Jumper,  Notre  Dame  U 
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•  Normalized  OPD 


OpDN0m  (/an/m)  = 


OPD 


RMS 


((/VAl  ) M  D) 


•  Elevation  angle  does  not  affect 
forward  looking  angles 

•  Large  increase  in  the  normalized 
OPDrms  at  large  backward  off- 
center  angles 

•  Centerline  is  a  “relatively  clean” 
region  of  the  flow  to  look  through 
at  backward  looking  angles 


Eric  Jumper,  Notre  Dame  U 


OPDNorm  (jjm/m) 
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Beam  Control  for 
Optical  Phased  Arrays 


Advantages  of  Phased  Arrays: 

•  Smaller  &  lighter 

-  An  array  of  independent  telescopes  takes  up  less  space  and  is 
lighter  than  a  monolith  and  beam  director 

-  Master  oscillator/power  amplifier  configuration  with  fiber  amplifiers 
takes  up  less  space  and  is  lighter  than  a  large  laser 

•  Conformal 


-  Telescope  array  is  conformal  to  the  surface 

-  Conformal  steering  elements  point  the  telescopes  individually 
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First  Target  Based  Phasing  Algorithm  that  Accommodates  Speckle 

The  4  Unobservable  Modes  Contain  Array  Tilt  and  Exponential  Class  Piston  Estimator 
Piston  Associated  with  3  Families  of  Aberrations  Accommodates  Unobservable  Modes 

Speckle 

Atmosphere 

Telescope 

0  5  10  15  20  25  30 

G-S  Iteration 


Speckle  Imaging  Synthetic  Aperture  with  Twice  the 
Resolution  of  the  Array  for  Aim  Point  Maintenance 


Stair  Mode  Imager  with  Matched  Filter  to  Correct  Array  Tilt 


Measurement  Discrepancy  Scales  with  Boresight  Error 

(Addresses  Multiple  Bean  Overlap  at  Target) 


Far  Field  Images  with  and  without  Stair  Mode 


Matched  Filter  Accurately  Predicts  Array  Tilt 
(For  SNR  =  1 ,  Strehl  =  0.9) 


RMS  1-Axis  Pointing  Error  (A7d) 


Glenn  Tyler,  tOSC 
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Speckle-Metric-Based 
Coherent  Beam  Combining 


•Propagation  of  an  HEL  beam  through  turbulence  leads  to  a  highly 
distorted  intensity  footprint  at  the  target  surface  (target  hit  spot). 

•The  hit-spot  beam  scattering  off  the  extended  target  results  in  strong 
speckle  modulation  at  the  HEL  beam  director  aperture. 


Experimental  results 

Hit-spot  with  speckle- 
metric  AO  OFF 


Jv 


Hit-spot  with  speckle- 
metric  AO  ON 


Fiber-array-based  HEL 
beam  director 


Target-return  speckle 
field  intensity 


Fiber-array1 


.based  HEL  b«ms 


Distorted 
target  hit 
spot 


/ 

/  i 

/  ^  m  I 

**  *  • 


Target-return 
speckle  field  phase 


The  target-return  speckle- 
field  has  highly  non-uniform 
intensity  and  phase  that  is 
composed  of  both  target-  and 
turbulence-induced 
components  with  a  large 
number  of  branch  points. 


Desired 
compensated 
target  hit  spot 


Svetlana  Lachinova,  Optonicus 
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Quantum  Key  Distribution  (QKD),  based  on  single-photon  quantum  states 

■Rapid  on-demand  generation  of  encryption  keys  among  air,  space,  and  ground 
communication  nodes 

■.Promises  provable  security  with  robustness  to  attacks  from  quantum  computing 
■  Challenge  is  overcoming  low  key  rates  in  free-space  transmission 


Approach: 

■  Utilize  spatial  modes  of  the  optical  field  to  achieve  high  keying  rates. 

■  Challenge  is  to  develop  techniques  for  encoding  and  decoding 
photons  with  spatial  modes  of  the  optical  field 


Accomplishment:  Validated  mode  decoding  with  transmission  volume  holograms 


hologram 
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in  ■  i 


phase 


m  =-2  -1 


incident  mode  number,  m 
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v  Summary  and  Goals 


•Observing  and  Identifying  Space  Objects 

•Improved  Imaging  of  Space  Objects 

•Improvement  of  imaging  capabilities  at  MSSS  and  SOR 
•Information  without  Imaging 

•Making  Space  Survillance  effective  without  large  telescopes 
•Predicting  the  Location  of  Space  Objects 

•Determination  and  prediction  of  many  orbits  at  high  accuracy 

•Remote  Sensing  in  Extreme  Conditions 

•Propagation  Through  Deep  Optical  Turbulence 

•Low  elevation  and  long  path  AO,  Active  Imaging,  and  Laser 
Propagation 
•Beam  Control 

•Good  beam  quality  and  imaging  at  higher  power 
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